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Abstract 



In this paper we develop a theory of mesoscopic fluctuations in disordered 
thin superconducting films in a parallel magnetic field. At zero temperature 
and sufficiently strong magnetic field the system undergoes a phase transition 
into a state characterized by a superfluid density, which is random in sign. 
Consequently, in this region, random supercurrents are spontaneously created 
in the ground state of the system, and it belongs to the same universality class 
as the two dimensional XY spin glass with a random sign of the exchange 
interaction. 
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Recent experiments on thin superconducting films in parallel magnetic field M have rekin- 
dled interest in this field. If the thickness of the films is small enough, the orbital effect of 
the magnetic field can be neglected and the suppression of superconductivity in the film is 
due to the Zeeman effect ' 2_4 l It has been observed that the resistance of such films at 
low temperatures and high enough magnetic fields exhibits very slow relaxation in timet 1 ]. 
This behavior is characteristic for spin and superconducting glasses. Below we discuss a 
possibility that mesoscopic fluctuations of superconducting parameters in disordered films 
account for such a behavior. Usually, in the limit where the electron elastic mean free path I 
exceeds the Fermi wave length h/kF, mesoscopic fluctuations of various physical parameters 
of superconductors are smaller than their averages ' 5 ~ 8 ' . Thus, they hardly affect macro- 
scopic observable quantities. However, there are situations where mesoscopic fluctuations 
determine macroscopic properties of superconducting samples. One example is a supercon- 
ductor in a magnetic field close to the upper critical field H c2 , where the magnetic field 
dependence of the superconducting critical temperature is determined by the mesoscopic 
fluctuations' 9 ' . In this paper we consider the case, where the magnetic field is parallel to the 
thin superconducting film and the main contribution to the suppression of superconductivity 
by the magnetic filed is due to Zeeman splitting of electron spin energy levels. We will show 
that at low temperatures T and high enough magnetic fields H, parallel to the film, the 
system exhibits a transition into a state where local superfluid density N s (f) (which is the 
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ratio between the supercurrent density J s and the superfluid velocity V s ) has random sign. 
In this case the system belongs to the same universality class as the two-dimensional XY 
spin glass model. The idea that the superfluid density can be negative has a long history 
[5-7,10-13] _ However, in the absence of magnetic fields and at zero temperature in slightly 
disordered superconductors (£o 3> I h/kp) the variance of the superfluid density, aver- 
aged over the superconducting coherence length £ = \Jd/Ao, turns out to be much smaller 
than its average! 5 " 8 ] ((5N S ) 2 } = G- 2 ({N S }) 2 < ((AQ) 2 , where 5N S = N s - (N s ), G is the 
dimensionless conductance of the normal metal film, in units of e 2 /h. Here D = vfI/3 is 
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the diffusion coefficient, A is the value of the order parameter at T, H — 0; vf is the Fermi 
velocity, and the brackets () denote averaging over realizations of random potential. This 
means, that as long as kpl h, the regions where the superfluid density is negative are 
rare and do not contribute significantly to macroscopic properties of superconductors. The 
situation in the presence of a magnetic field parallel to the film is different, because the 
average superfluid density decays with H faster than its variance. Hence, at high enough 
magnetic field the amplitude of the mesoscopic fluctuations of N s (r) becomes larger than 
the average, and the respective probabilities of having positive and negative signs of N s (r) 
are of the same order even at kpl 1. 

A theory of magnetic field induced phase transition, which does not take into account 
mesoscopic fluctuations predicts I 2 ' 14 ' 15 ] that at low temperatures the superconductor-normal 
metal transition is of first or second order depending on whether the parameter A r so is larger 
or smaller than unity, respectively. Here r so is the spin-orbit relaxation time. 

Let us start with the case where A r so « 1. At T = and within an approxi- 
mation which neglects mesoscopic effects, the value of the critical magnetic field H® is 
the result of the competition between the average superconducting condensation energy 
density (E c ) ~ ^oA 2 , and the polarization energy of the electron gas in the magnetic 
field. Here vq is the average density of states in the metal on the Fermi surface. The 
average spin polarization energy density of nonsuper conducting electron gas is of order 
of (E p (0)) ~ v (fi B H) 2 . Its relative change in the superconducting state is of order of 
(E p (0)) - (E p (A)) ~ §7rA r so (£ p (0)) < (E p (0))^ 18 l As a result we get an expression 
for the critical magnetic field H° = H cc (Aqt so )~i 3> H cc . Here H cc = A //i_B is the 
Chandrasekar-Clogston critical magnetic field of the superconductor-normal metal transi- 
tion for A t so — > oo and Hb is the Bohr magneton. 

Now let us consider the mesoscopic fluctuations of the quantities, discussed above, in a 
volume whose size is of the order of the coherence length £ - The amplitude of mesoscopic 
fluctuations of the polarization energy is of order of t 19 ' yj ((5E P ) 2 ) ~ G^ 1 (A t so ) (E p (0)) , 
while its A-dependent part is of the order of 
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y/((6E p (A(H))-6E p (A = 0))*) ( A(H)\ 2 

Ur • (1) 



Here (A(H)) = A Q ^J (H® — H)/ is the average superconducting order parameter. Since 
both the polarization energy and the condensation energy are fluctuating quantities, A(r) 
should also be spatially fluctuating. Let us consider a domain of size Ld £o where the 
value of A(r) differs from its bulk value by a factor of order of unity. An estimate for the 
energy of such a domain consists of three terms, namely 
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M(A) 1 Co , r H° c -H g 

v A(H)*dL* D {Cl GL^ + ° 2 ~ Hf~ + Cz LV (2) 

where d is the thickness of the film and Ci,C 2 ,C 3 are factors of order of unity. The first 

term in Eq.2 corresponds to the A-dependence of mesoscopic fluctuations of polarization 

energy and has a random sign. When estimating this term we have taken into account that 

domains of size £ make independent random sign contributions into Eq.2. The second and 

third term are the average condensation energy and surface (gradient) energy of the domain, 



respectively. It follows from Eq.2 that if L D ~ £(H) = £oy H®/\H — H®\ that there is an 
interval of magnetic fields near the critical one H® — H ~ H®/G 2 , where the first term is 
larger than the second and the third ones. It means that, in this case the spatial distribution 
of A(r) is highly inhomogeneous and the amplitude of the spatial fluctuations of A(r) is of 
order of its average, while the characteristic size of the domains is of order of L D ~ £(H = 
H®(1 — l/G 2 )) ~ £ G. Superfluid density in this region has a random sign as well. To see 
this, one should consider states with finite superfluid velocity V s = (V% + 2e/cA)/m, where 
x{t) is the phase of the order parameter, A(r) is the vector potential of a magnetic field 
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perpendicular to the film and m is the electron mass. If V s (r) is of the order of the critical 
velocity, all three terms in Eq.2 are modified by factors of order of unity when compared with 
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the case V s = 0. The second and the third term in Eq.2 decrease with V s , while the first term 
is changed in a random direction. This means that at high enough magnetic fields, states 
with nonvanishing value of V s (r) have lower energy than the states with V s = 0, and that 
the system is unstable with respect to the creation of supercurrents of random directions. 



In this estimate we neglected the energy of the magnetic field associated with V s (r). Since 
at each point of the system the possible energy gain associated with finite value of V s (r) 
is independent of the direction of V s , the ground state of the system is highly degenerate 
and belongs to the same universality class as XY spin glass with random sign of exchange 
interaction. 

It is important to mention that even in the case of small magnetic fields in the presence 
of spin orbit scattering the time reversal symmetry is broken and the electron wave functions 
are complex. These currents flowing in the random directions exist even in normal metals. 
By evaluating the diagrams shown in Fig. 2a, we derive the correlation function of the current 
density in normal metals(|r — r'\ ^> h/k F ), 

{.m-W)) « Siij^TTsoipLsHySir- P). (3) 

Here r = 1/vf is the elastic mean free time and i,j are coordinate indexes. It is important 
to note, however, that for a given configuration of the scattering potential and at a given 
value of the external field the spatial distribution of J(r) is a unique function. This implies 
that the currents described by Eq.3 do not exhibit features which can be associated with 
superconducting glass states. 

Below, we will be interested in supercurrents much larger than those described by Eq.3. 
Such currents are spontaneously created at strong enough magnetic fields as a result of 
the instability associated with the random sign of superfluid density. Consider the Gorkov 
equation for A(f)' 20 ', 

A(f) = g J dr'Kif, P; H, A(r), {A(f)})A(r'), (4) 

where K(r, P; H, A(f), {A(r)» = kT £ e Gf(r, P; H, A(f), {A(f)})^GT e (f, P; H, A(f), 0) 
a* a ; G"P(r, f 1 ; H, A(f), {A(r)}) is the exact one particle Matsubara Green function, a, (3, v, \x 
are spin indexes, a x a ^ is the Pauli matrix and e = (2n+l)iikT is the Matsubara frequency, g is 
the dimensionless interaction constant. Both A(r) and K(f, r') in Eq.4 are random functions 
of realizations of scattering potential in the sample. Averaging Eq.4 over realizations of the 
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random potential and using the approximation (A(r)K(f, r')) = (A(r; H))(K(r, r', H)) we 
get the above mentioned expression for H®. In the case of strong magnetic fields, when 
A(r; H) <C A , we can expand Eq.4 in terms of A(r). Since A(r) varies slowly over distances 
of the order of £ 0) while (K(r,r f )) decays exponentially for \ f — r'\ 3> £o, we can also make 
the gradient expansion of Eq.4. As a result we get from Eq.4 

( 1 2p - H° — H\ r -, -, -, A 3 (r) 

(j2& 2 (V - ~ A ? + A ^ + 7 ^V>',^)A(r')dr' = °(-^), (5) 

where 5K°(r,r') = K°(f,P) - (K°(f,?)) and K°(r,r') = K(f,r',{A(r) = 0}). The dif- 
ference between Eq.5 and the conventional Ginsburg-Landau equation is the third term in 
Eq.5 which accounts for mesoscopic fluctuations of the kernel K°(f,r'). It is precisely this 
term, which at high magnetic fields leads to the random sign of superfluid density. 

Employing the perturbation theory with respect to 5K°(f, r') we get from Eq.5 an expres- 
sion for the correlation function C(fi — r 2 ) = (5A(fl)5A(f* 2 )) of the mesoscopic fluctuations 
of the superconducting order parameter SA(f; H) = A(r; H) — (A(H)) 

c«-v~§\ \n-^«my, (6) 

G [ exp(-2 v / 3|n - f 2 |/£(#)), |rl - f 2 | » £(H). 

In order to derive Eq.6 we had to calculate the correlation function (SK°(f{, rtySK ^, r^)) 
using the diagrams shown in Fig.2b. It follows from Eq.6 that {(5A) 2 ) in the two-dimensional 
case is almost independent of H, but (A(if)) decreases with H. As a result, perturbation 
theory holds as long as (A(H))/A = y '(H° - H) / ' H° c > G' 1 . 

Using the expression for the supercurrent expanded in terms of A(r, H) <C A we have 
for the correlation function of the nonlocal superfluid density N s (f,r'), 

2 2 

JJjf) = — ( df'N s {r- r')V s {?) = — [ dr'{N s {H)5{r- ?) + 5N s (r,r'))V s (r'), 
m J m J 

(iV^ (W" (n ~ ^ " " 2) ^r[F (n " ( ^ ~ " 2) ( } 

which is valid as long as (A(H))/A > G' 1 and £(if) > |n - f 2 | > £ - At 
|r*i — r 2 | ^> £(H), the correlation function in Eq.7 becomes exponentially small. Here 



N S (H) = N°(A(H)) 2 /A 2 , JVj = N(l/£ ) 2 is the average superfluid density at H = and N 
is the electron concentration in the metal. The first term of the correlation function in Eq.7 
is connected to the fluctuations of the order parameter A(r) as shown in Fig.2c. The second 
term corresponds to the fluctuations of the Green functions G e (f, f) shown in Fig. 2d. 

Therefore if the magnetic field is close to the critical one, i.e. \H — H®\/H® ~ G~ 2 , then 
the amplitude of fluctuations of the superfluid density averaged over the size £(H) becomes 
of order of its the average 5N S ~ (N s ), which means that the local value of the superfluid 
density, averaged over the size £o> becomes of random sign. Hence the system is unstable 
with respect to spontaneous creation of supercurrents. 

If \H — H®\/H® <C G~ 2 , one can neglect the second term in brackets in Eq.5. Rescaling 
f ~ xG£ 0: A(r) ~ A /Gf(r/G£ ), yields a dimensionless stochastic equation for f(x) 



where {Sk(x,x')) = and the correlation function {5k(x\, x' 1 )5k(x 2 , x' 2 )) = (S(xi — x[) + 
G~ 2 /{G~ 4 + {xi — x' 1 ) 4 })5(a ? 1 — x 2 )5(x[ — x' 2 ) is given by diagrams shown in Fig.2b. It 
follows from Eq.8 that the amplitude of fluctuation of the modulus of the order parameter 
SA(r) ~ (A(H)) ~ A /G is of order of its average. The characteristic spatial scale of the 
fluctuations of 5A(r) is of order of Ljj. The sign of second term in Eq.8 fluctuates randomly 
which corresponds to the random sign of the superfluid density. The spontaneously created 
supercurrents in this case have random directions, their typical amplitude is of order of 
Jg ~ N®h/G 3 !;o and their characteristic scale of spatial correlations is also of order of Lp. 
The current described by Eq.3 is negligible compared with when I <^ G^ 1 ^/ Dt so . 

The fact that the sign of N s is random is especially clear in the case of large magnetic 
field, when H — H® ^> H®G~ 2 . In this case, A(r) is nonzero only due to existence of the rare 
regions, where 5k(x,x') is much larger than the typical value. Thus, the spatial dependence 
of the modulus of the order parameter has the form of superconducting domains embedded 
in a normal metal. These regions are connected via the Josephson effect. We can calculate 
the average critical current of the junctions and its variance as functions of the distance 




(8) 
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between the droplets L : 

They decay with L exponentially and as a power law respectively. As a result, the amplitude 
of the fluctuations turns out to be larger than the average, which means that J c has a random 
sign. 

It is well known' 21 ! that at T = the long range order of the ground state of the two- 
dimensional XY model is destroyed by an arbitrary small concentration of antiferromagnetic 
bonds. As we have mentioned above in the case H <C regions, where N s (f) < 0, exist with 
small but finite probability. In this case, however, the properties of superconducting system 
are different from the XY model because the supercurrents spontaneously created in these 
regions are screened by the Meisner effect. Thus at H, T = superconducting films should 
exhibit the conventional long range order. This implies that there is a critical magnetic field 
H S g < H® where at T = the system has a phase transition from superconducting to the 
superconducting glass states (if° — Hsg ~ H®G~ 2 ). The interval of magnetic fields where 
the system is in the superconducting glass state is indicated by shaded region in Fig.l. 

Let us now consider the case of weak spin-orbit scattering limit A r so ^> 1. In this case 
the spin magnetization in the superconducting phase is zero. Correspondingly, the conven- 
tional theory based on the equation for average order parameter leads to the conclusion that 
the superconductor-normal metal transition is of first order with the critical magnetic field 
HcJ 3 ' 4 ^. However, the fluctuations of both magnetization energy of the normal metal and 
the condensation energy of the superconductor phase should lead to a nonuniform state, 
qualitatively similar to the case A r so <C 1. The theory of this phenomenon at A r so 1 
is, however, more difficult. In this case a domains of normal phase within a bulk supercon- 
ductor (or a superconducting domain in normal metal) has the surface energy of order of 
dL d^ou AI, where Lp is the domain size. This energy is much larger than the typical energy 
associated with mesoscopic fluctuations in Eq.2, dL^QUQ^sH^G^ 1 . Thus the probability 
of the occurrence of such domains is small as long as G > 1. We would like to stress though, 
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that qualitatively the case A r so ^> 1 is not different from the case A t so <C 1 for in both 
cases the superconducting glass solutions survive at T = and H > H®. 

The question whether or not the quantum fluctuations of the phase of the order parameter 
destroy the superconducting glass state at T = and large H is still open^ 22-24 ! . At finite 
temperatures T > 0, strictly speaking, the system considered above doesn't posses a phase 
rigidity because of Meisner screening effect ' 2 ^. On the other hand the two dimensional XY 
model with random sign of exchange interaction is known not to exhibit a phase transition 
between the paramagnetic and the spin-glass phases' 2 ^ , again implying the absence of long 
range order. We acknowledge useful discussions with B.Altshuler and S.Kivelson. This 
work was supported by the Division of Material Sciences, U.S. National Science Foundation 
under Contract No.DMR-9625370 and the US-Israeli Binational Science Foundation grant 
no. 94-00243. 
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FIG. 1. Qualitative picture of the magnetic field dependence of A(H) at zero temperature 
when Aot s o <C 1. The shaded region corresponds to the superconducting glass phase. A* ~ Ao/G, 
5H C ~ tf c °/G 2 . 
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FIG. 2. a) Diagrams representing the current correlation function Eq.3. b) Diagrams rep- 
resenting the correlation function (SK°(f{, fl)SK°(f2, rjj)). c)d) Diagrams representing the cor- 
relation function of supercurrent densities (J s (r)J s (r')}. Solid lines correspond to electron Green 
functions in metal and dashed lines correspond to elastic scatterings of a random potential and 
black dots represent the correlation function (SA(fi)SA(f2)) given by Eq. 6. 
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